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Due to efficient scaling with electron number N, density functional theory (DFT) is 
widely used for studies of large molecules and solids. Restriction of an exact mean-field 
theory to local potential functions has recently been questioned. This review summarizes 
motivation for extending current DFT to include nonlocal one-electron potentials, and 
proposes methodology for implementation of the theory. The theoretical model, orbital 
functional theory (OFT), is shown to be exact in principle for the general N-electron 
problem. In practice it must depend on a parametrized correlation energy functional. 
Functionals are proposed suitable for short-range Coulomb-cusp correlation and for long- 
range polarization response correlation. A linearized variational cellular method (LVCM) 
is proposed as a common formalism for molecules and solids. Implementation of nonlocal 
potentials is reduced to independent calculations for each inequivalent atomic cell. 



1. Introduction 

Hohenberg and Kohn 1 (HK) proved for nondegenerate ground states that total 
electronic density determines the external potential acting on an interacting N- 
electron system. The constrained search procedure of Levy 2 extends this proof to 
general variational forms of N-electron wave functions, easily extended to spin- 
dependent densities. Spin indices and summations are assumed here, but omitted 
from notation. HK theory defines a universal ground state density functional F s = 
E — V, which reduces for noninteracting electrons to the kinetic energy functional 
T s of Kohn and Sham 3 (KS). Implementation of the implied density- functional 
theory (DFT) followed only after introduction of the KS orbital mean-field model. 
In KS theory, the density function is p(r) = X^ n i/°i( r )' a sum of orbital densities 
Pi = 4>*4>i for the model state, with occupation numbers rij. The many review 
articles and monographs on DFT are exemplified by Parr and Yang 4 and by Dreizler 
and Gross 5 . 
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2. Motivation for an exact orbital theory 

2.1. Can DFT with local potentials be an exact theory? 

The relationship of DFT to many-body theory remains obscure, requiring quite 
different methodology for strongly-correlated systems. The KS construction is not 
easily related to strong-correlation theory, in which a Hubbard model is often in- 
troduced that has no relationship to the electronic density. There is no system- 
atic way to improve model functionals. Nonphysical self- interaction produces in- 
correct ground states for magnetically ordered systems. Nonbonding interactions in 
molecules and optical potentials in electron scattering are not correctly described. 

Numerous problems or paradoxes in DFT are associated with the assumption 
that exact mean- field theory requires only local one-electron potential functions 6 . 
It is important to note that the mathematical issue here is not the existence of 
HK functionals as such, but rather the existence of the uniquely denned density 
functional derivatives required to determine variational Euler-Lagrange equations 7 . 
Quantitative tests show that an optimized effective potential (OEP local exchange) 
does not reproduce ground-state energies and densities in the unrestricted Hartree- 
Fock (UHF) variational model for atoms 8 ' 9 . Because the UHF energy functional 
is uniquely denned, and the accurately computable UHF ground-state density is 
both noninteracting and interacting v-representable 5 , this textbook example of con- 
strained variational theory is a rigorous consistency test for DFT 7 . The exchange- 
only limit of DFT linear response theory 10 is inconsistent with the time-dependent 
Hartree-Fock theory of Dirac 11 , due to failure of locality 12 . A postulated local ki- 
netic energy would imply an exact Thomas-Fermi theory (TFT), in conflict with 
KS DFT and with the exclusion principle for more than two electrons 13 . 

In an orbital model, such as KS theory, the exclusion principle is imposed 
by independent variation and normalization of each occupied orbital of a model 
state 7 . This logic produces an extended TFT 14 when applied to orbital densities. 
This orbital-density theory is operationally equivalent to orbital functional the- 
ory (OFT), and reduces to KS DFT in the local density approximation (LDA) 15 . 
There is no conflict with rigorous theory 16 . Confined to normalized total densities, 
functional analysis implies neither Thomas- Fermi nor KS equations 15 . 

2.2. OFT with nonlocal potentials can be an exact theory 

For nondegenerate states, the KS Ansatz, p = ^2 i niPi, expresses p as a function 
of occupied orbitals of a Slater determinant model state $. With this Ansatz, 
any well-defined density functional is also an orbital functional. If F s [p] = E — 
V is parametrized as an explicit functional of p, it defines an orbital functional 
F[{(/>i}] for all functions in the orbital Hilbert space. Schrodinger variational theory 
is expressed in terms of orbital functionals with functional differentials of the form 7 

6F = Y^nij*T{6ti-^+cc}. 
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This defines an orbital functional derivative of the general form — J 7 <!>%■> 

where T is a linear operator. These functional derivatives determine orbital Euler- 
Lagrange (OEL) equations Tfa = {ti — v(r)}(f>i. This is an exact theory if the true 
N-electron F = E — V can be expressed as an orbital functional 17 . 

OFT is deduced from formally exact N-electron theory and nonlocal exchange- 
correlation potentials are constructed. Quite explicitly, the exchange "potential" is 
the Fock exchange operator, and the correlation " potential" is determined by orbital 
functional derivatives of the linked-cluster expansion 7 . The implied OEL equations 
are a realization of "exact" KS equations. OFT implies a systematic refinement 
of the LDA based on implementation of nonlocal potentials. The theory is not 
density-dependent, bypassing questions of v-representability and of the convergence 
of density-gradient expansions. It is not restricted in principle to ground states, to 
weak correlation, or to short-range interactions 17 . It systematizes many alterna- 
tives currently being pursued in molecular DFT: the orbital-dependent exchange 
originally proposed by KS 3 , and the OEP approximation 7,8 , sometimes character- 
ized as "exact KS" 18 . OFT nonlocal potentials are derived from functionals defined 
by antisymmetric wave functions. In principle, this eliminates self-interaction and 
the need for any ad hoc correction 19 . In practice, computational methodology de- 
veloped for the latter may be very useful in the context of the orbital-indexed 
potentials considered here. The present analysis has important practical implica- 
tions for computational methodology relevant to the electronic structure of large 
molecules and solids. With a new theoretical basis embedded in N-electron theory, 
extension of current work on exact exchange in DFT to include electronic corre- 
lation requires timely development of efficient computational methods appropriate 
to nonlocal potentials 20 . 



3. Orbital functional theory 

3.1. An exact mean-field theory exists 

Specify any relationship * — > $ that determines a model state $ for any N-electron 
state ^, where $ is a Slater determinant for any nondegenerate state. This postu- 
late implies the KS Ansatz in a variational theory of the model density. Imposing 
normalization (3>|\I r ) = (3>|$) = 1, which causes no formal problems for a finite 
system, (H - = implies that E = E Q + E c , where E Q = ($|i7|$), and 

E c = E — Eq = (&\H \ty — $) defines the correlation energy. Since Eg is an explicit 
orbital functional, an exact mean-field theory exists if E c can be derived from many- 
body theory as an orbital functional. The linked-cluster expansion shows that this 
is possible. Exact OEL equations are defined formally by orbital functional deriva- 
tives of the linked-cluster expansion. An implicit closed form of the exact functional 
E c has been derived 21 . 
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3.2. Kinetic energy and exact exchange 
u 

2 

T = ^2ni(i\t\i); V = ^2m(i\v\i) 



Using i = — tjV 2 and v = v(r), functional E = T + U + V is defined such that 



for orthonormal orbital functions. The Coulomb interaction u = l/r\2 defines or- 
bital functional U = Eh + E x , where 

E h = \'Yl,'riinj{ij\u\ij)\ E x = ~'Y^nin j {ij\u\ji). 

i,3 i,j 



3.3. Variational equations for the orbitals 

Orbital Euler-Lagrange (OEL) equations are derived using the orbital functional 
derivatives 

- , SU . , SV 

—ju = u fc> = v \ r m- I 1 ) 



niS<j)* ' n z 8cj>* ' niScj)* 

Here u = Vh(r) + v x , where Vh is the Coulomb potential, and v x is the Fock ex- 
change operator. An orbital functional E c [{4>i}] defines ^§pr = v c 4>i- It will be 
assumed that some parametrized E c is defined such that v c is Hermitian. A La- 
grange functional consistent with independent orbital normalization is defined by 
subtracting m[{i\i) — l]ei from the energy functional. Defining F = E — V, the 
OEL equations are 

F<j>i = {d - v}<j>i, i = l,---,N, (2) 

using the orbital functional derivative n S ^, = T4>i = {i+ u + v c }4>i. Here i and u 
are explicitly nonlocal (linear operators), and there is no implication that v x + v c 
is equivalent to a multiplicative local potential. 



4. OFT as an exact variational model 

4.1. Correlation energy as an orbital functional 

Unsymmetrical normalization ($|\I> — $) = defines an orthogonal projection, 
= $ and = # - where V = and Q = I - V. The N-electron 

Schrodinger equation can be partitioned by projection. In the limit — > 0+, 

= -{Q(H -E -E c - i^Q]- 1 ^. 

Hence, for 77 — * 0+, 

e c = ($|#|* - *) = (*|ir|e*) 

= -(*\H[Q(H -Eo-E c - i^Q^H^). 

Because model state $ determines Q, this exact expression for E c is an implicit 
orbital functional 21 . 
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4.2. Gateaux functional derivatives 

Suppose that an orbital functional F is denned such that = Vf<f>i and Of 

is Hermitian in the orbital Hilbert space. With no additional assumptions, the 
functional differential 



SF = Y / n i J rf 3 r{^I(r)«/^(r) + cc}-^n 4 J d 3 r 



determines a density functional derivative. This derivative is an orbital-indexed 
local potential function 

5F = (3) 



pi 

defined throughout the orbital Hilbert space 20 . 

For unconstrained orbital variations about any stationary state, the functional 
differential of F = E — V is determined by Eqs.(2), the OEL equations. In detail 22 , 

SF = E n * / d 3 r{5cf>*(r)FMr) + cc} = ]T n t J d 3 r{e, - v(r)}5 Pi {r) 

i i 

determines the density functional derivative 

SF r\ n 

— — =v fi {r) =ei-v(r). 

nidpi 

Hence for any stationary state, Vfi(r) depends on the orbital subshcll index unless 
all orbital energies are equal 22 . Such dependence on a "direction" in the Banach 
space of densities defines a Gateaux functional derivative, the generalization to 
functional analysis of an analytic partial derivative 23 . The Gateaux derivative ex- 
pressed by an indexed local potential Vfi(r) is in general singular at nodes of the 
orbital function <pi(r). A Frechet derivative 23 (multiplicative local potential) can 
exist only if vji is independent of its orbital index at all coordinate points, which 
requires elimination or exact cancellation of the orbital nodes. This condition is 
imposed by the construction of v xc in the LDA, but cannot be assumed to be a 
general consequence of variational theory. 

In order to generate Euler-Lagrangc equations with a normalization constraint, 
a variational functional of total density must be defined for arbitrary infinitesimal 
variations about normalized densities 7 . The HK functional F s is defined only for 
normalized ground states. No Euler-Lagrange equation is implied unless this def- 
inition can be extended to include infinitesimal neighborhoods of ground states. 
If a total (Frechet) functional derivative were to exist, it would justify such an 
extension, implying an exact Thomas- Fermi theory 15 . If the functional differential 

SF = [ d 3 r S -fs P 
J op 

is a unique functional of the total density variation dp, it defines a Frechet functional 
derivative 

^=v f (r) (4) 
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as a multiplicative local potential function. Comparison of Eqs.(4) and (3) shows 
that the total functional differential is not unique unless the corresponding Gateaux 
derivatives are all equal. For variations about a ground state, this requires all orbital 
energies to be equal, violating the exclusion principle for any compact system with 
more than one electron of each spin. Hence exact TFT is not implied for more 
than two electrons. In contrast, the Gateaux derivatives are well-defined, and imply 
modified Thomas-Fermi equations with orbital indices, consistent with the exclusion 
principle, and with the general form of KS equations . 

Eqs.(3) and (4) cannot be reconciled for the HK functional F s unless all or- 
bital energy eigenvalues are equal 13,15 . This difference is a consequence of different 
normalization constraints in DFT or OFT and in TFT. That independent normal- 
ization of orbital densities is physically correct, specifically enforcing the exclusion 
principle in Hartree-Fock theory and KS DFT, is obvious for the example of an atom 
with two noninteracting electrons of parallel spin. The ground state is ls2s 3 S'. The 
density constraints are: 

DFT : J p u = 1, J p 2s = 1 

TFT: j(p ls +p 2 s) = 2. 
This TFT constraint allows the nonphysical solution 

TFT : J p ls = 2, J p 2s = 0, 
which violates the exclusion principle. 



5. How to do it 

The technical problem of generalizing KS DFT to an exact theory reduces to 
methodology for indexed local potentials. For atoms or local atomic cells, this just 
adds an indexed correlation potential w C i(r) = to standard Hartree-Fock 

methodology for the Fock exchange potential 24 . Off-diagonal Lagrange multipliers 
can be used to constrain orbital orthogonality during self-consistent iterations. For 
closed-shell systems, these Lagrange multipliers must vanish on convergence. 



5.1. Energy-linearized variational cellular method 

Basis-set expansion becomes impractical for large molecules and solids. The alter- 
native methodology of multiple scattering theory, used in many variants for the 
electronic structure of condensed matter 25 , can be formulated as a variational cel- 
lular method (VCM) 26 ' 7 . In full-potential MST, local basis functions are computed 
numerically in each atomic cell of a space-filling Voronoi lattice. These functions 
are truncated at local cell boundaries, and linear combinations are matched varia- 
tionally across cell interfaces 7 . The logical structure and inherent efficiency of the 
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VCM emulate a tight-binding model. It can be linearized in energy and adapted 
to a general variational method for solids and molecules 27 . Using numerical basis 
functions in local cells validates semirelativistic calculations for heavy atoms, and 
validates frozen-core approximations. These features are exploited in condensed- 
matter methodology. Truncating basis functions at cell boundaries removes inter- 
cell overlap and multicellular Coulomb integrals from the methodology. Local cells 
interact through boundary matching and through extended potential functions that 
act as an external field in each cell. As in quantum field theory, the essential physics 
is that of quasiparticles interacting through effective mean fields. 

5.2. The local nonlocal (LNL) model 

Intracell calculations including an indexed correlation potential are just like 
Hartree-Fock for atoms, not a practical problem. Indexed local potentials are con- 
sistent with MST, which simply uses a different orbital basis in each local cell. 
The residual practical problem is the treatment of long-range tails of the indexed 
potentials. The truncated basis removes all intercell overlap contributions to these 
potentials. Electronic correlation arises from two distinct causes. The short-range 
effect of the Coulomb cusp is significant for intracellular correlation, while the long- 
range effect of polarization is expected to be the principal intercellular correlation 
effect. In the case of metals, long-range correlation is required for dynamical ex- 
change screening, needed to eliminate nonphysical behavior of exact exchange at 
the Fermi surface. The long-range asymptotic form of the indexed potentials for 
exchange and correlation is electrostatic, modified by polarization response. The 
significant success of DFT and of model theories suggests that, not far outside a 
source charge density, the asymptotic potentials in general reduce to nonindexed 
collective fields. If this is valid, it would justify a "local nonlocal" (LNL) model, 
in which true indexed potentials are used only within local cells, while asymptotic 
forms of the extracellular fields are combined as an effective mean field in each such 
cell. As in the Ewald expansion, this implies significant cancellation or screening of 
fields. This may justify neglecting the multipole part of the effective field in each 
cell during self-consistent iterations, including it as a first-order perturbation after 
convergence. 

5.3. Coulomb-cusp correlation 

Short-range correlation is dominated by the Coulomb pole u = I/V12, requiring 
a f + \r\2 cusp in the wave function. The Colle-Salvetti (CS) Ansatz 28 imposes 
this cusp condition, multiplying model state <& by a symmetrical factor IL<j [1 — 
£(rj,rj)], where, in relative coordinates q = r« — rj and r = i(r^ + r^), £(r,q) = 
exp(— [3 2 q 2 )[l — r(r)(l + \q)] 29, ■ Correlation energy can be expressed exactly as a 
sum of pair-correlation energies, E c = (<f>|_ff|5' — $) = J2i<j n i n j e ij- This can be 
parametrized, in analogy to CS, by e?- = —{ij\u£ij{q)\ij), where u = u(l — V12), 
P12 is the 2-electron exchange operator, and £ij(q) = exp(— /3?-g 2 )[l — 7ij(l + \q)}- 
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The normalization condition (<j>|\& — $) = determines 

_ fa"|exp(-^g 2 )|ij) 
113 ~ (y|exp(-/^.g2)(i+ig)|y)' 

such that = for each pair. The free parameter (iij can be chosen 

to minimize the Bcthc-Goldstonc (independent electron pair model) energy for 
each orbital pair. The nonlocal correlation potential is defined by v c <j)i = n .f^* = 
~ n j{j\^ij\j)4'i- Antisymmetry is built in. Terms j = i vanish, so there is no 
self-interaction. The indexed local correlation potential is 

v ci (r) = -^^^rijijlu&jlj^iir). 

j 

5.4. Multipole response correlation 

If the mapping ^ — ► $ is determined by the Brucckner-Brenig condition 29 ' 30 , 
= 0, so that $ is the reference state of maximum weight in this re- 
moves all single virtual excitations <&" from \P. Then the leading terms of Q\& are 
double virtual excitations <f> a „ b with coefficients c" h . Higher-order terms do not in- 
teract with <E> through the Coulomb interaction. Thus E c = ($\H\Qty) is a sum of 
pair correlation energies 7 , e£- = J2 a <b(^ ~ n a)(l — nb)(ij\u\ab)cfj. The expression 
(a|t) c |i) = (a\5E c /ni5<f>*) can be evaluated, treating the coefficients cfj as constants 
in the current cycle of an iterative loop 21 . These matrix elements 

(a\v c \i) = ^rij ^(1 - n c )(l - n b )(aj\u\cb)(cb\c\ij) 

j c<b 

- ^2 uktij ^2(1 - n b )(kj\u\ib)(ab\c\kj), 

k<j b 

agree with the leading self-energy diagrams in many-body perturbation theory. This 
determines the kernel of the linear operator v c . 

As an illustrative example, this formula has been simplified and applied to the 
leading long-range term in an iterative expansion of the multipole polarization po- 
tential implied by bound-free correlation between a scattered electron and a polar- 
izable target 21 . This analysis introduces a polarization pseudostate <j> Pj , computed 
as the normalized first-order response of an occupied orbital <j>j of the target sys- 
tem to an external multipole field. For multipole index A > 0, transition matrix 
elements (pj\u\j) vary as r A for small r and as l/r A+1 for r > ro. The implied 
multipole polarization potential is quadratic in these transition elements and varies 
as r 2A and l/r 2A+2 , respectively, in these limits. Hence, unlike earlier models from 
purely asymptotic theory, this potential vanishes as r — > 0. When applied for A = 1 
to variational calculations of low-energy e-He scattering, this polarization poten- 
tial produced scattering phase shifts in good agreement with earlier variational 
calculations that incorporated dipole polarization response. The indicated general 
methodology is to compute first-order multipole pseudostates within each atomic 
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cell, and then use them to construct the multipolc polarization potentials that 
contribute to the mean external field in all other cells. 

6. Conclusions 

This brief review is intended to motivate and to outline methodology that should 
be implemented in order to get beyond current practical limitations of density- 
functional theory. It is a call to the condensed-matter community to move beyond 
DFT by extending current methodology for nonlocal potentials. It is a call to the 
theoretical chemistry community, concerned with the electronic structure of large 
molecules, to set aside CI methods and restricted DFT in favor of an efficient, 
practicable, and ultimately exact ab initio mean-field theory. 
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